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ABSTRACT 



Aims. Ultra compact HII (UCHII) regions are indicators of high-mass star formation sites and are distributed mainly in the Galactic 
plane. They exhibit a broad band spectrum with significant emission between near-IR and radio wavelengths. We intend to investigate 
the possible contribution of the forthcoming ESA Planck mission to the science of UCHII regions by evaluating the possibility of 
detecting UCHIIs that are bright in the radio regime. 

Methods. We performed new 7 mm observations of a sample of UCHII regions. The observations were designed to acquire high- 
frequency radio spectra. For each source in our sample, the free-free radio spectrum has been modeled. Along with far-IR measure- 
ments, our spectra allow us to estimate the flux densities of the sources in the millimeter and sub-millimeter bands. We extrapolated 
and summed the ionized-gas (free-free radio emission) and dust (thermal emission) contributions in the afore mentioned wavelength 
ranges. The possibility of Planck detecting the selected sources can be assessed by comparing the estimated flux densities to the ex- 
pected sensitivity in each Planck channel. To obtain a realistic estimation of the noise produced by the Galactic emission, the Planck 
sky model software package was used. 

Results. For each target source, from our new 7 mm data and other radio measurements from the literature, important physical 
parameters such as electron density and their spatial distribution, source geometry and emission measure were derived. We conclude 
that, in the case of the present sample, located close to the Galactic center, Planck will have a very low detection rate. In contrast, 
assuming that our sample is representative of the whole UCHII-region population, we derive a very high probability of detecting this 
kind of source with Planck if located instead close to the anticenter. From the analysis of the ionized-gas properties, we suggest that 
the selected sample could also be contaminated by other kinds of Galactic objects. 

Key words. H II region - Radio continuum: ISM - Stars: formation - ISM: dust, extinction - Stars: circumstellar matter 



1. Introduction 

The Planck mission of the European Space Agency (ESA) 
will operate over a wide range of frequencies, from 30 to 900 
GHz (from 1 cm to about 300 //m). The frequency coverage 
of the Planck channels was selected to correspond to the so- 
called cosmic microwave background (CMB) window, where 
temperature fluctuations due to Galactic emission are minimal. 
Although Planck was designed for cosmological studies, its 
multi-frequency nearly-full-sky maps will have a profound im- 
pact on fundamental physics as well as Galactic and extragalac- 
tic astrophysics and seve ral Galactic projects are indeed in cluded 
in its core program (see lThe Planck Collaboration] d2005l) for de- 
tails about the scientific case for the Planck mission). 

In the science of compact Galactic objects, Planck will be 
sensitive to the millimeter emission from dusty envelopes of 
young and old stars and to the thermal radio emission from the 
ionized gas surrounding many classes of stellar objects. 

HII regions of size about 0. 1 pc, emission measure higher 
than 10 7 pc cm 4 , electron density higher than 10 4 cm" 3 , and 
ionized mass of around 10~ 2 M Q are classified as ultra compact 
HII regions (UCHIIs) and are sites of high-mass star formation. 
In these regions, newly born stars are still embedded in their na- 
tal molecular clouds. The ambient dust thermally reradiates the 
entire luminosity of the central star in the infrared (IR), mak- 
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ing UCHIIs among the brightest IR sources in the Galaxy. The 
radiation fields of the massive hot stars inside UCHIIs (OB spec- 
tral type) ionize the circumstellar material. The signature of ion- 
ization is the typical free-free flux-density distribution at radio 
wavelengths (S v °c v"). Therefore, we expect the Galactic popu- 
lation of UCHIIs to be a strong foreground contamination source 
to the Planck experiment. 

As part of the preparation of a Planck pre-launch catalog, we 
started an observing program to measure the 7 mm flux density 
of samples of Galactic sources with t he 32 m INAF-IRA Noto 
Radiotelescope. In a previous paper (Um ana et al.ll2008al) . we 
published the our flux density measurements of planetary nebu- 
lae (PNe). 

In this paper, we present new 7 mm single-dish observations 
of UCHIIs. Main goal was to obtain reliable estimates of the 
flux densities expected in Planck channels by building and mod- 
eling the spectral energy distribution (SED) of our targets. The 
observing frequency (43 GHz) is within the cosmic microwave 
background (CMB) window and represents one of the channels 
of the forthcoming Planck mission. Therefore, in this frequency 
band, we obtain a direct measurement of the expected flux den- 
sity with minimal frequency extrapolation. 

As added value, the present work provides a sizeable dataset 
of 7 mm measurements of UCHII regions. These data provide 
strong constraints on the observed SEDs in the spectral region 
where free-free emission dominates. 
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Table 1. Observations. 
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Reported as a multi source in the 6 cm catalog but associated with the same MSX source. The 7 mm observations were carried out by pointing 
at the MSX coordinates. 



The SED modeling in the radio-millimeter spectral range is 
a crucial step for the study of the physics of ionized envelopes 
around UCHII regions. The 7 mm measurements presented in 
this paper provide the possibility of constraining parameters 
such as the electron density and the radial profile, the source 
geometry, and the emission measure. 



2. Sample selection 

Our sa mple was selected from the catalog of iGiveon et al.l 
(2005a), where a cross-correlatio n between the M idcourse Space 
Experiment (MSX6C) catalog dEgan et al.ll2003l) an d the Multi- 
Array Galactic Plane Imaging Survey (MAGPIS) dWhite et alj 
2005) was performed. MAGPIS was carried out at 6 cm and 
covered Galactic regions with longitude 350° < / < 42° and lat- 
itude \b\ < 0.4°. A comparison of the two catalogs (radio and in- 
frared) allowed the authors to find a sample of 687 radio sources 
with highly-reliable infrared counterparts. Only 2 28 of 687 tar- 
gets w ere detected in the new MAGPIS at 20 cm (Helfan d et al.1 
120061) . 

All sources with both 6 and 20 cm data have red MSX col- 
ors and flat or inverted spectral indices, as is typical of thermal 
radio sources. Th e average spectral ind ex between 20 and 6 cm 
is equal to 0.28. IGiveon et al.l d2005at) classified these sources 
as compact HII regions on the basis of their radio and infrared 
properties. For most sources, this was their first classification. 
Starting from this subsample of 228, we first inspected all the 
corresponding MAGPIS 6 cm images. 



To minimize the source-structure effect on the flux den- 
sity measurements, we selected only sources compact enough 
to be considered point-like at the resolution of the Noto 
Radiotelescope at 7 mm (HPBW about 1') and rejected those 
located in high-confusion regions. In our sample, we included 
single sources of 6 cm angular size < 20" and sources that, al- 
though mapped as multiple at 6 cm, showed separations < 10" 
and were associated with only one infrared counterpart. Finally, 
from the resulting sample we selected only sources with inte- 
grated flux densities higher than 100 mJy at 6 cm and declination 
higher than -15°. In the conservative hypothesis of an optically 
thin nebula at 6 cm, a cut-off at both 100 mJy and 6 > -15° 
guarantees 7 mm detectability and observability with the Noto 
Radiotelescope. The expected 7 mm flux density is a lower limit 
in the case of an optically thick nebula at 6 cm. All the above 
selection criteria reduced our sample to a total of 51 UCHIIs. 

The names and positions of the selected targets are reported 
in Cols. 1, 2, and 3 and 6, 7, and 8 of Table Q] When observing 
multiple sources associated with the same infrared counterpart, 
we pointed the telescope at the MSX coordinates. 



3. Observations 

The observations reported in this paper were carried out at dif- 
ferent epochs between January 2007 and July 2008, using the 
32 m Radiotelescope in Noto (Italy). The telescope has an active 
surface system that accounts for the gravitational deformation of 
the primary mirr or, providing high qualit y observations at high 
frequencies (see lOrfei et aD d2002l 120041) for details of the ac- 
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tive surface system mounted on the Noto radiotelescope). The 
7 mm receiver used in our observations is a cooled superhetero- 
dyne receiver of typical zenith system temperature (r sys ) about 
100 K. The observations were performed with a 400 MHz instan- 
taneous band centered on the sky frequency v s k y = 42950 MHz, 
the beam size being 0.9'. The gain ranges from 0.04 to 0.08 K/Jy, 
depending on the air mass as a function of the elevation. Our tar- 
gets were observed using the On The Fly (OTF) scan technique, 
which consists of driving the beam of the telescope across the 
sources in the RA direction. The typical scan duration was ap- 
proximately 20 s, which is short enough to remain close to the 
white-noise regime of the radiometer. The scan length was 3', 
corresponding to about three telescope beams. To achieve a good 
signal-to-noise ratio, each source was observed several times, for 
a total integration time of 30 minutes. Multiple OTF scans were 
then added together. 

Daily gain curves were obtained and the flux scale fixed us- 
ing NGC 7027 as a primary calibrator and 3C286 and DR21 as 
secondary calibrators. The adopted flux density for NGC 7027 
is 5.07 Jy . This value was obtained by correcting that reported in 
lOtt et al.l d 19941) for t he observed evoluti on of the spectral flux 
density in NGC 7027 (IZiilstra et alj|2008h . 

All measurements were corrected for atmospheric extinc- 
tion. The correction factor was calculated by performing a series 
of tipping scans at elevations between 20° and 80°, covering a 
range of sec(z) from 1.01 to 3. Tipping scans were performed 
several times (once per hour) during each observing run. The 
opacity correction was applied to each flux measurement after 
interpolating the results of contiguous tipping scans. The sky- 
opacity measurement averaged over all of the observing sessions 
was about 0.14. The opacity fluctuation ranges from about 7% 
(measured during observing sessions with very stable sky condi- 
tions) to about 30% (observing sessions characterized by worse 
weather). 

The telescope pointing was checked before each observing 
session by measuring its accuracy for a large sample of pointing 
calibrators: SiO maser sources and bright sources at 7 mm. The 
estimated pointing accuracy is about 10". 

Taking into account the telescope pointing accuracy, the er- 
ror in the flux density of the primary gain calibrator, the errors in 
the opacity measurements, and the gain-curve error, we estimate 
the overall calibration accuracy to be around 10-15%. 

We detected all targets with a 3cr confidence threshold, with 
a detection rate of 100%. The results of these observations are 
reported in TableQ] where the measured 7 mm flux densities (S 7) 
are listed with their errors (Cols. 4, 5 and 9, 10). 

4. The SEDs 

4.1. The dust contribution 

To determine the dust contribution for our sample of UCHIIs, we 
firstly searched both the literature and online catalogs for data in 
the range from 8 /mi to about 1 mm, where the dust contribu- 
tion i s relevant. As a co nsequence of the selection criteria, MSX 
data (Ega n et aT1l2003l) are availabl e for all of the 51 so urces 
of the sample. We found IRAS data dHelou & Walkerll 1988b for 
45 source s, Scuba data for 12 sources at 450 jxm and for 16 at 
850 / mi (IWilliams et al.ll2004t IWu et al.ll200l iThompson et"ai1 
2006), and Calthech Submillimeter Observa tory (CSO) d ata for 
4 sources at 350 /mi dHunter et all 120001: iMotte et all 120031) . 
For 8 sources, we found measurements at 1.2 mm carried out 
with the IRAM 30 m Radiotelesc ope dAltenhoff et all 119941: 
Beut her et all |2002|) and SIMBA (Faun dez et all 120041) . Ten 
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Fig. 1. SED template of the dust contribution. The data collected 
from the far-IR to the mm range have been normalized to 1000 
Jy at 100 /mi. The squares represent the mean values of the nor- 
malized measurements and standard deviations are also shown. 
Two dust components are evident: the hot component determines 
the contribution at shorter wavelengths (dotted line), while the 
cold component dominates in the submillimeter spectral range 
(dashed-line). 



sourc es were also obser ved at 1.3 mm at the IRTF on Mauna 
Kea (IChini et al.lll986allbl) . at the SEST dLaunhardt & Henningl 
119971) . and with the MAMB O bolometer array installed on the 
IRAM 30 m Radiotelescope (IMotte et alj|2003l) . To summarize, 
we have enough measurements to compile SEDs for 34 sources 
in our sample. From the SEDs, we derive an SED template by 
normalizing the 100 yum high-quality IRAS data to 1000 Jy and 
then averaging together all the individual SEDs. 

The SED template obtained can be modeled by using two 
dust components at different temperatures. We assume that each 
of the two dust components emits as a blackbody modified by 
the frequency-dependent dust opacity. This gives S v °c v^B Y (Td), 
where B V {T^) is the Planck function for a dust temperature of 
7d and where B is the emissivity index. The index B strongly 
depends on the mineralogical composition of th e grains and on 
t heir p hysical shapes. Following the results in Chini & Kriigel 
(1986) about modeling dust emission from UCHII regions, we 
find that for the cold dust component a good fit is achieved for a 
dust temperature of 30 K and an emissivity index B — 2. The hot 
component has been fitt ed for a dust temper ature of 120 K and 
by assuming that B — 1 (Faundez et al. 2004). 

In Fig.Q] we show the averaged data and overplot their best- 
fit solution. The deviation of the point at 8 /mi could be related to 
the contribution of either a very hot dust component or the stellar 
component, whereas the deviation of the two measurements at 
1.2 and 1.3 mm can be explained as evidence of an ionized gas 
contribution. 

For 17 sources, we do not have enough measurements at 
A > 100 /mi to compile their individual SEDs. To estimate their 
SED contribution from cold dust, we assume that all sources in 
our sample have the same dust properties. Our SED template is 
therefore a good assessment of their emission curves (Fig.[TJl and 
we can obtain an individual SED for each of them by just scal- 
ing the template to their MSX 21 /mi measurement. In the SED 
template, the flux density at 21 /mi is found to be a factor of 0.04 
of its value at 100 /mi. 
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Table 2. Sources parameters. 
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0.11 


0.11 


31.0702 + 0.0498 


0.34 


0.19 


16.9440 - 0.0736 


0.23 


6.17 


26.6085 -0.2122 


-0.09 


0.26 


31.2803 + 0.0615 


0.01 


0.21 


18.4608-0.0034 


0.02 


3.39 


27.4938 + 0.1893 


0.3] 


0.03 


30.6671 -0.3316 


0.13 


4.77 


18.7104 + 0.0001* 


0.45 


> 12.80 


27.9777 + 0.0780 


0.27 


0.11 


31.2430-0.1108 


0.16 


1.12 


18.3029-0.3910 


-0.13 


0.11 


28.2461 + 0.0134 


0.24 


0.10 


32.1514 + 0.1317 


-0.14 


0.06 


18.7084-0.1263 


-0.03 


3.76 


28.2007 - 0.0494* 


0.58 


> 2.21 


31.3948-0.2585 


0.53 


0.13 


19.7403 + 0.2799 


0.07 


0.03 


28.6096 + 0.0170 


0.21 


3.25 


32.4727 + 0.2041 


0.34 


3.12 


18.8330-0.3004 


0.41 


0.50 


28.6520 + 0.0271 


0.09 


0.61 


32.7977 + 0.1903* 


0.31 


> 3.41 


21.8732 + 0.0076 


-0.03 


0.13 


28.2875 - 0.3639 


-0.06 


1.89 


32.2718 -0.2260 


-0.15 


0.25 


21.3853-0.2537* 


0.32 


>2.49 


29.2113 -0.0690 


0.08 


0.73 


33.8104-0.1869 


0.32 


0.03 


22.7253 -0.0106 


0.28 


0.07 


29.9564-0.0174 


0.09 


3.00 


35.4672 + 0.1381 


0.04 


0.76 


23.9548 + 0.1494 


-0.10 


0.12 


30.5345 + 0.0209 


0.08 


1.24 


35.5736 + 0.0678 


0.19 


0.26 


24.5076 - 0.2222* 


0.35 


> 0.76 


30.2340-0.1392 


-0.15 


3.16 


35.5789-0.0304* 


0.62 


> 8.96 


25.5196 + 0.2156 


0.04 


0.29 


30.8667 + 0.1141 


0.34 


5.38 


38.5497 + 0.1636 


0.03 


0.03 


25.2666-0.1584 


-0.09 


10.71 


31.4134 + 0.3092 


-0.13 


0.56 


38.8750 + 0.3088 


0.08 


2.25 


25.3977-0.1406 


0.00 


0.60 


30.6877 - 0.0729 


0.06 


1.72 


39.8821 -0.3457 


-0.07 


1.75 


25.3820-0.1813 


0.29 


0.80 


30.7206 - 0.0826 


0.11 


3.76 


41.7410 + 0.0972 


-0.09 


0.08 



For these sources the simulated radio spectra indicate that the emission is optically thick up to 7 mm, the derived EM has to be considered as 
a lower limit. 



4.2. The ionized-gas contribution 

We estimate the contribution from ionized gas to the SEDs of the 
selected sources by com bining our 7 mm me asurements with the 
flux densities at 6 cm in lGiveon et a l. (2005a) and those at 20 cm, 
the observations at the lower frequency being available for each 
source in our sample from the new MAGPIS catalog reported by 
iGiveon et al. I (l2005bl) . 

In the present work, we need to assess the contribution 
from free-free emission in an angular element of size com- 
parable to the Noto Radiotelescope beam (HPBW about 1'). 
For those sources with multiple peaks in their 6 cm and 
20 cm high-resolution observations, we use as an estimation 
of their 6 cm and 20 cm flux densities the sum of the in- 
tegrated fluxes of each peak. Literature data at radio wave- 
lengths shor ter than 6 cm we r e also collected for 1 8 targets in 
our sa mple dWink et al.lll982h IWood et alJ I 1988alfbl: TKurtz et all 
1994 iGehrz et alJ Il995t iFev et all 119951: IWalsh et all Il998t 
Martin-Hernandez et al]|2003l:lLongmore et al.ll2007l). 



The spectral index a between 6 cm (IGiveon et alJ20 05a) and 
7 mm (this paper) was calculated for each source. The values 
obtained are reported in Table [2] The free-free emission from 
the ionized gas was derived by modeling all the available radio 

measurements for e very source. 

As described by lHoare et alJ dl99 lh . a typical UCHII region 
can be modeled using spherical symmetry. The source can be 
modeled as a central cavity surrounded by a shell of ionized gas, 
the newly born high-mass star being located inside the cavity 
center. The source geometry is characterized by the linear di- 
ameter and the ratio of the internal to external radii (ri nt /r ext ). 
The shell of ionized material is located in the inner region of 
the dust cloud responsible for the infrared emission. In this sim- 
ple spherical-geometry model, the spatial density structure of the 

ionized gas is described by the power law n e {r) oc r~ y . 

Ad opting the model software developed by lUmana et alJ 
(2008b), the free-free radio spectrum was simulated for each 
source by assuming a fixed value of the source linear diameter. 
In the following, we consider two cases of the source linear di- 



ameter (equal to 2xr ext ): 0.1 pc and 0.05 pc. The free-free model 
spectrum was calculated for different values of the ratio nntAext, 
the index of the power law that describes the ionized-gas spatial 
distribution, and the corresponding average electron density. 

The distance to the source was treated as a free parameter in 
obtaining the best-fit SED to both the observations and simula- 
tions. The source angular size derived from the assumed linear 
diameter and distan ce needs to be comp arable to the angular size 
measured at 6 cm (IGiveon 



s to be compa 
et al.ll2005a). 



The new 7 mm measurements allow us to characterize well 
the free-free radio spectrum of each source in our sample. The 
distance-independent parameters are, in fact, univocally defined. 
The values of y and rj nt /r ext for the two sizes do not show any 
significant difference, as shown in Fig. [2] Using our model we 
found that the index y covers wide range of model solutions, 
many sources being reproduced well by assuming homogeneous 
electron-density radial distribution, while others have distribu- 
tions that are highly inhomogeneous. In all cases, the geometry 
of most of the sources could be described by the ratio r mt /r ext 
below 0.2. The derived electron density n e is on the order of 10 4 
cm 4 . 



Following Terzian & Dickevl d 1973b . it is possible to derive 
the mean emission measure (EM) if the observed radio optically 
thin flux density (S 7, in mJy) and the corresponding source an- 
gular size (6, in arcsec) are known, 



EM = 



5.3 x 10 5 S 7 

e 2 



[pc cm 6 ] 



(1) 



To derive the mean emission measure from our 7 mm measure- 
ments, we assume that the source ang ular size at 7 mm e quals 
the angular size at 6 cm measured bv lGiveon et alJ d2005a). The 
derived EM is reported in Table [2] For sources with simulated 
spectra that are optically thick up to 7 mm, the derived EM is 
considered to be a lower limit. 

Note that the spectral index between 6 cm and 7 mm by it- 
self is not able to characterize the emission at 7 mm. In many 
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Fig. 2. Comparison between the derived model parameters for 
the two analyzed source diameters (0.1 pc and 0.05 pc). Left 
panel: index of the power law that describes the density gradient 
inside the source. Right panel: ratio of internal to external source 
radii. 



cases, sources with positive spectral index are characterized by 
optically thin emission at 7 mm. 

Our modeling also allows us to estimate the emission mea- 
sure. In Fig. [3] the emission measure obtained from the model fit 
is plotted as a function of the EM calculated from Eq. Q] 

The model used to simulate the radio emission from UCHII 
regions assumes the highly simplified hypothesis of spherical ge- 
ometry, whereas radio images at high angular resolution show 
that these sources are chara c terized by various morphologies 
dWood & Churchwelll [T989at iKurtzet al.lll994t iDePreeetaT] 
2005), i.e., shell-like, bipolar, cometary, irregular, and spheri- 
cal. In spite of the possible difference between the real and the 
modeled source geometry, the clear correlation between the val- 
ues of EM deduced from the model and those calculated from 
the observations indicates that our model parameters allow us to 
reproduce well the overall properties of the targets. 

Our calculations exhibit a large spread in EM (Fig. [3). HII 
regions characterized by EM > 10 7 pc cm -6 are defined as ultra 
compact. Sources with EM values at the edges of the range do 
not show any significant difference in their dust emission com- 
pared to the other sources of the sample. No relation was ob- 
served between infrared colors and EM. This physical param- 
eter is also randomly distributed as a function of the Galactic 
longitude. 

The infrared colors of the selected sample follow a typical 
trend for UCHII regions, while in the radio-millimeter regime 
not all sources have typical UCHII characteristics, as illustrated 
by their EM values. Our study in the radio range demonstrates 
that a sample of UCHIIs selected only on the basis of their in- 
frared colors may be contaminated by other classes of Galactic 
sources. Planetary nebulae or high-mass evolved stars could 
have infrared colors similar to those of UCHII, but different ion- 
ized gas properties. 



5. Estimated detectability with Planck 

Planck is equipped with a Low Frequency Instrument (LFI), op- 
erating at 30, 44, and 70 GHz and a High Frequency Instrument 
(HFI) operating at 100, 143, 217, 353, 545, and 857 GHz. To 
evaluate the fraction of UCHIIs in our sample that is detectable 
with Planck, we need to compare the modeled flux densities of 
the sources with the expected sensitivity in each observational 
band. 

We derive the modeled flux densities at each Planck channel 
from the source SEDs by calculating the dust and gas emission 
as described in Sects. |4.1| and |4~2l and then summing both con- 



B 10.00 r 




S 0.10 r 



.01 0.10 1.00 10.00 100.00 

EM/10' [pc cm - *] (7 mm obs.) 

Fig. 3. Comparison of the Emission Measure deduced from the 
simulation with those calculated from the 7 mm measurements. 
Arrows refer to those sources optically thick ut to 7 mm. 



tributions. In Fig. [4] four of the derived SEDs are shown as an 
example; the shaded area represents the spectral range covered 
by the Planck satellite. The dust contribution was obtained by 
extrapolating the modified-blackbody fit of the cold dust compo- 
nent (dashed lines in Fig. |4j, while the free-free contribution was 
obtained by modeling the data points as described in Sect. 14.21 
and extending the derived synthetic spectra to the submillimeter 

range (dotted line in Fig. [Hi . 

Following lUmana et al.l (120061) . the expected total sensitivity 
is assumed to be the sum in quadrature of all the sources of con- 
fusion noise, namely (i) the nominal instrumental Planck sen- 
sitivity per resol ution element, (ii) th e extragalactic-foreground 
confusion noise dToffolatti et al.ll 1998b . (iii) the CMB anisotropy 
dBennett et al.ll2003a ). and (iv) the Galactic -foreground confu- 
sion noise. 



As pointed out bv lWood & Churchwelildl989bl) . the Galactic 
distribution of candidate massive-star formation sites is mainly 
confined to the Galactic plane, and our sample consists indeed 
of UCHIIs located in the Galactic plane, where we also expect 
the Galactic emission to be the major source of confusion noise 
in all Planck channels because of the emission of the dust and 
gas located in the Galactic disk. 

The Galactic noise reported by Uma na et al. I d2006h was de- 
rived for regions of the sky outside of the Galactic plane, there- 
fore we cannot use this value for UCHII regions. For a real- 
istic estimation of the noise related to the Galactic emission, 
we simulate the sky emission in the range between 30 and 900 
GHz with the Planck sky model (PSM) software package (ver- 
sion 1.6.3). PSM is a set of programs and ancillary data devel- 
oped by the Planck ESA Working Group 2 that can simulate 
the sky emission at Planck frequencies and angular resolution, 
and reproduces the thermal free-free, synchrotron, and the dust 



components, i.e., thermal dust an d spinning dust. 



component is modeled following lDickinson et al. 



31 The f ree-free 
d2003|)andthe 

spatial structure estimated from the dust-extinction corrected H a 
emission. The synchrotron component is calculated by assum- 
ing a perfect power-law with different values of the spectral in- 
dex in different sky directions, as derived from the 408 MHz 
(lHaslam et al.lfT982l) and WMAP 23 GHz dBennett et al.ll2003bh 
maps. Th e thermal emission fro m interstellar dust is assessed 
following Fin kbeiner et al.l (fl999). The prediction of the rotating 



1 See Delabro uille et alj (2009) and lLeach etal] d2008h for applica- 
tions and capabilities of the PSM package. 
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Fig. 4. Examples of the source SEDs, different symbols refer to 
different components. Diamonds: hot dust component; squares: 
cold dust component; triangles: ionized gas. The dust component 
is deduced from the template SED, the free-free emission of the 
ionized gas derived from its model. The spectral range covered 
by Planck is represented by the shaded area. 



grain emission is obtained b y multiplying the co lumn densities 
of the cold neutral medium dSchlegel et alJ ll998) with the emis- 
sivities calculated using the spinning dust model developed by 
iDraine & Lazarianl (119981) . 

For the two Planck channels of lower frequency (30 and 44 
GHz), the Galactic emission is dominated by the free-free com- 
ponent, whereas for the third channel (70 GHz) the free-free and 
the thermal dust components are comparable; for the other chan- 
nels the emission of the thermal dust dominates. The other emis- 
sion mechanisms provide negligible large-scale contributions. 

To obtain an idea of the trend of the Galactic emission fluc- 
tuations as a function of the Galactic longitude, the standard 
deviation of the Galactic emission for each PSM channel map 
was calculated in boxes of size 3 x V (6 V is the Planck beam at 
the considered frequency), centered on Galactic latitude 0° and 
spaced every 5° in longitude. We selected the size of the region in 
which the standard deviation of the Galactic emission was to be 
calculated as a compromise between the necessity both to have 
a large number of pixels in the region and to remain well inside 
the plane of the Galaxy. We show the RMS of the Galactic emis- 
sion for each channel map at different values of Galactic longi- 
tude for the two Planck instruments in Fig. [5] The Galactic emis- 
sion fluctuation strongly depends on the longitude. The first and 
fourth Galactic quadrants (center) are characterized by higher 
noise with respect to the second and third quadrants (anticenter), 
as an obvious consequence of the larger amount of dust in the 
direction of the center of the Galaxy. 

The correctness of the PSM predictions was tested by com- 
paring them to the observed all-sky WMAP temperature map for 
Stokes I at 41 GHz (Q band). The WMAP data were retrieved 
from the NASA Legacy Archive for Microwave Background 
Data Analysis (LAMBDA), while a sky temperature map at the 
same frequency and angular resolution was constructed with 
PSM. We found very good agreement between them, allowing 
us to use our results on the Galactic emission fluctuation with 
confidence. 

In Fig. [6] we show the percentage of objects from our sample 
that could be detected above the 3<x threshold for each channel of 
the instruments onboard Planck. Unfortunately, all the UCHIIs 



Fig. 5. Expected Galactic emission fluctuation for \b\ < 0.5° at 
different values of Galactic longitude for the two Planck instru- 
ments. 



in our sample are located in the first and fourth quadrants and 
because of the high Galactic emission fluctuation a very small 
fraction of them should be detected by Planck. 

However, UCHIIs have been detected in any Galactic direc- 
tion (IWood & Churchwelllll989bl) . Therefore, assuming that our 
sample is representative of the entire class of UCHIIs, we can 
also calculate the detectability of those located in the region of 
the Galactic anticenter. Figure [6] shows how UCHIIs in the di- 
rection of the anticenter should be easily seen by Planck. 

Unfortunately, the high-mass star-formation rate is expected 
to strongly decrease f ar away from the center of the Galaxy. 
iBronfman et al.1 d2000t) reported that the OB star-formation sites 
in the anticenter direction reduced to about 10% those detected 
in directions close to the center of the Galaxy. The expected 
increased Planck capability to detect UCHIIs in the anticenter 
direction has to be compared with the decreasing star forma- 
tion rate. The true expected number of detectable UCHIIs must 
be weighted by the radial distribution of star formation in the 
Galaxy. 

The UCHII sample analyzed in this paper repr esents only a 
subsam ple of the 228 UCHII regions reported by iGiveon et alj 
(2005a), which are distributed in a region of the Galactic disk 
that cover about 52° of Galactic longitud e. Taking into account 
the estimation of IBronfman et al.l d2000l) . the expected UCHIIs 
that Planck could detect for a similar range of Galactic longi- 
tudes in the anticenter direction could be tens. Note that the 
Planck channels that have the higher capability of detecting 
UCHIIs are the channels at the lower wavelengths (Fig. [6), 
where the dust emissi on contribution is d ominant. However, 
the sample selected by IGiveon et al. I (I2005al) has a bias towards 
sources above the threshold detectability at the radio wave- 
lengths. Therefore, it is reasonable to expect that the estimation 
reported here is only a lower limit to the true number of UCHII 
regions detectable by the Planck satellite, as there may be many 
sources with very low radio emission levels, for actual sensitivi- 
ties, but with detectable IR dust emission. We emphasize that this 
conclusion is supported by many detections of UCHIIs in the re- 
gion of the Galactic anticenter obtained with the ball oon exper- 
iment Archeops at mm and sub-mm wavelengths dDesert et al.l 
120081) . 

In many cases UCHIIs have been observed on the borders of 
larger and fainter HII regions. In higher wavelength channels, the 
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Fig. 6. Expected Planck capability of detecting the sample of 
UCHII analyzed here. Dashed line: true source Galactic posi- 
tion. Dot dashed line: results for the same population but located 
in the anticenter region. 

Planck beam may be large enough to detect emission from these 
associated normal HII regions. Considering their contribution, 
the Planck capability to detect UCHIIs here reported could be 
an underestimate. 

We conclude that we expect the Planck mission to provide an 
important contribution to studies of UCHIIs, as it will provide 
the scientific community with measurements in a still largely- 
unexplored spectral region. 

6. Summary 

We have presented 7 mm (43 GHz) observations of a sample 
of radio-bright ultra-compact HII regions carried out with the 
INAF-IRA Noto Radiotelescope. The results of this survey pro- 
vide a sizeable database of new millimeter measurements that 
allow us to characterize the high-frequency free-free contribu- 
tion from the ionized-gas components of the selected targets. 
The free-free spectrum of each source were derived by modeling 
all the available radio and millimeter observations and the new 
7 mm data. 

We used our 7 mm measurements along with infrared and 
sub-millimeter data to compile SEDs from the radio to the far-IR 
and to estimate the expected flux densities in the spectral region 
between the cm and the sub-mm ranges. 

The study of properties of Galactic UCHIIs is one of 
the secondary science objective s of the Planck core program 
dThe Planck Collaboration|[2005h therefore, it appears to be im- 
portant to test the capability of Planck to detect UCHII regions. 
We, thus, estimated the fluctuation of the Galactic emission as 
a function of the Galactic longitude by using the Planck Sky 
Model software package. In this paper, we note that we report 
for the first time a quantitative estimation of the emission fluctu- 
ation of the Galactic disk as seen by Planck. As expected, the 
Galactic emission fluctuation is higher near the center of the 
Galaxy. When comparing the modeled radio/sub-mm flux den- 
sities to the expected total sensitivity of the Planck satellite, we 
evaluate that few sources are likely to be detected by Planck in 
the direction of the Galactic center, where the sources of our 
sample are located. If we assume that the sample analyzed in the 
present paper is representative of the whole Galactic UCHII pop- 
ulation, we infer a higher rate of detectability by Planck in the 



direction of the Galactic anticenter, where the Galactic emission 
fluctuation is significantly lower. 

Therefore, even though designed for cosmological studies, 
Planck should also contribute to the UCHII science. The rich 
source database provided by Planck will allow follow-up obser- 
vations aimed at the investigating morphological details with fu- 
ture instruments working in the same spectral region, such as 
ALMA. 

The source parameter emission measures were derived both 
directly from the 7 mm measurements and, independently, from 
a model of the ionized-gas continuum emission. The derived EM 
values for each sources are clearly correlated, supporting the 
goodness of the model fit. The analysis of this parameter leads 
us to conclude that our sample could be contaminated by other 
types of Galactic sources, although all of the selected targets can 
be classified as UCHII regions because of their infrared colors. 
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